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Abstract: In this research work the flow behavior was analyzed, which affect the power coefficient as well as torque 
coefficient of a helical Savonius rotor is investigated by means of commercial code CFD. Conventional three bladed 
Savonius rotors have high coefficient of static torque at certain rotor angles and negative coefficient of static torque at 
certain rotor angle. In order to increase the efficiency of all rotor angles the system a helical Savonius rotor with a twist of 
certain proposed degree is introduced. A three-dimensional model of 90 D twist three blades helical Savonius rotor at 

0 D , 60 D , 120 c and 150 : rotor angle has been constructed by using the software gambit of the Fluent 6.3 package. The 
contours of static pressure and velocity magnitude around the rotor blades area at horizontal iso -plane is obtained from the 
CFD simulation. High performance was obtained at advanced bucket in upstream airflow at 60 D rotor angle and maximum 

positive static pressure obtained at 150 c rotor angle, which affect the positive coefficient of static torque. 
Keywords: Helical Savonius rotor, rotor angle, CFD analysis, velocity contour, and static pressure contour. 



I. Introduction 

The Savonius rotor is a vertical axis wind turbine (VAWT) originally invented by a Finnish engineer Sigurd J. 
Savonius in 1922 [1]. N. Fulisawa observed the relation between rotor performance and flow characteristics field around a 
Savonius rotor by using a laser light sheet and a CCD camera with a shutter speed of 0.002s [2]. Nobuyuki Fujisawa 
measured the flow fields in an around Savonius rotors at various overlap ratios with particle imaging velocimetry technology 
and comparison analysis is conducted between a vertical-axis spiral rotor and a conventional Savonius rotor by using CFD 
technology [3]. The flow field around a Savonius rotor, which is closely related to the torque and the power performance of 
the rotor, has been investigated by many researchers to clarify the power mechanism of the rotor. Most investigations have 
been experimental studies with flow visualization techniques, such as tufts [4, 5]. The flow physics of the rotor with overlap 
variation 12.37% to 25.87% were analysed with the help of velocity and pressure contours of the rotor. From this analysis 
the overlap of 19.87% was the optimum overlap condition at which pressure and velocity differences across the rotor were 
the highest for which maximum power extraction by the rotor would be possible at that overlap condition [6]. The 
combination of Darrieus and Savonius rotors thus exhibits high starting torque and high power coefficient. Keeping the 
above points they have made to study the flow physics of a combined three -bucket Darrieus and Savonius rotor using Fluent 
6.0. They have considered four overlap conditions between 10.87% to 25.87% for the Savonius part of the combined rotor 
and observed highest relative speed and highest pressure of the combined rotor near its top at 15.37% overlap condition [7]. 
M.A.Kamoji et al. [8] performed the static torque coefficients at all the rotor angles for helical Savonius rotors are positive. 
However, for conventional Savonius rotor, there are several rotor angles at which static torque coefficient is negative. 
Conventional Savonius rotor configurations with steady two-dimensional RANS CFD simulations were also conducted by 
Gupta et al. in which static torque predictions for a Savonius rotor with external flow Control demonstrated good agreement 
with the corresponding experimental wind tunnel results [9]. Gupta et al. [10] has also been observed that the aerodynamic 
performance and optimum overlap of the combined rotor are strongly controlled by the flow physics of the individual rotors. 
Katsuya I et al. [11] simulated numerically the flow over a two -bucket Savonius rotor by solving 2D finite volume equations 
using an upwind scheme for unstructured mesh. The result shows good matching with the experimental values of power 
coefficient, aerodynamic coefficients. M., Manzoor Hussain, Mehdi et al. [12] conducted the enhancement of efficiency by 
modifying the blades configuration from straight semicircular to a twisted semi circular one. Wind flow analysis was done 
over each configuration of the rotor with the blade twist ranging from 5 p to 60 D in a steps of 5 C . Maximum efficiency is 

obtained 33.85% at 45° twist angle, when compared to 25.6% without twist. A fully automatic optimization finally takes 

place, the commercial tool Gambit for geometry and grid generation (including quality check) and the industrial CFD code 
ANSYS-Fluent to compute the flow field around the Savonius rotor [13]. 

II. Description Of Physical Model 

In the present research work, the performance of a three -bucket helical Savonius rotor with semicircular profiles 
buckets at 90° twist angle as shown in Figure- 1 was analyzed. The operation of the rotor is simulated with the variation in 
rotor angle atO^ 1 , 60 D , 120^ and 150 c . The diameter of the rotor has a value of 24cm, while the value of the rotor height is 
40cm. The thickness of the buckets is 3 mm. 
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III. Details Of Computational Domain 

The three dimensional computational domain of three -bladed helical Savonius rotor along with the boundary 
conditions is shown in the Figure- 2. The computational domain resembles was 60cm x 60cm in cross -section. Wind speed 
was considered on the inlet face ABCD, outflow was considered on the outlet face EFGH, symmetry condition was 
considered for the side walls of the computational domain and blades are consider as rotating wall. The computational 
domain was adopted of the unstructured hexahedral mesh and the moving mesh techniques in Figure- 3. The computations 
were initially carried out with various levels of refinement of mesh which is done in gambit of the Fluent 6.3 package. The 
number of cells and nodes involved in meshing of the rotor with various levels are given in Table 1 . The grid independent 
limit computational mesh (Masson eta al., 1997) [14] was obtained by changing the resolution of mesh at all important areas 
of the rotor which is shown in the Figure- 4. 

3.1 Computational Domain Discretization Method: 

The continuity equation, Newton's second law, the conservation of momentum equation (N-S equation) and 
realizable k-s turbulence model are numerically solved by the commercial software FLUENT 6.3. Finite Volume Method is 
used for discretizing a differential equation and computes the static pressure and velocity fields of the flow on the rotor 
buckets. Applying the discretizing techniques the computational domain divided into a finite number of elementary control 
volumes. The sequential algorithm, Semi-Implicit Method for Pressure-Linked Equation (SIMPLE), was used for solving all 
the scalar variables with the k-s turbulence model. The finite volume discretization transforms the differential equations 
governing the flow in a linear system of algebraic equations that is iteratively solved. 

3.2 Mathematical Formulation: 

In the present study, the results were obtained by applying the Finite Volume Method to solve the conservation 

equations of the turbulent air flow on the helical Savonius rotor, which allows the calculation of the performance of the rotor. 

Among these conservation equations, there is the mass balance equation. The Eq. (1) with indicial notation represents the 

— -) 

mass balance, where p is the density, V is the velocity vector and value of 5„. is zero for steady-state flow. 

■£+V.( P V) = S m = 0 (1) 

The momentum equation must be solved along with the mass balance equation. The Eq. (2) represents the 
momentum equation, where p is the density, V is the velocity vector, p is the static pressure, and T is the stress tensor 

(described below), and p5 and F are the gravitational body force and external body forces (e.g., that arise from interaction 
with the dispersed phase), respectively. 

j: (pvy v. (pi? V) = - 17 P + v.T+pg +F (2) 

The stress tensor is given by 

f = u [(V.F + V.V 7 ) - z-V.Fl] (3) 

*Zi 

Where u is the molecular viscosity, I is the unit tensor, and the second term on the right hand side is the effect of volume 
dilation. 

3.3 Turbulence Model: 

The realizable k-s turbulence model provides superior performance for flows involving rotation, boundary layers 
under strong adverse pressure gradients, separation, and recirculation. The modeled transport equations for k and 8 in the 
realizable k-s model are: 




Where 



www.ijmer.com 



740 I Page 



International Journal of Modern Engineering Research (IJMER) 
www.ijmer.com Vol.3, Issue.2, March- April. 2013 pp-739-746 ISSN: 2249-6645 

?7 k. ■ 

C 1 = max [O.43.73T], r| = S.- and S = / 25, { S t ,- 

"T j ^* j 

In these equations, G k represents the generation of turbulence kinetic energy due to the mean velocity gradients. 
G b is the generation of turbulence kinetic energy due to buoyancy. represents the contribution of the fluctuating 
dilatation in compressible turbulence to the overall dissipation rate. Cn And C 1 _ are constants. CT k and cr_ are the turbulent 
Prandtl numbers for k and 8 , respectively S k and 5 _ are user-defined source terms [15]. 



IV. Figures and Tables 




Figure 3: Computational mesh of 90° twist helical Savonius rotor 
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Table 1: Details of refinement levels of the 3 -buckets helical Savonius rotor 



Refinement 
Level 


No. of Nodes 


No. of Elements 


1 


34226 


179994 


2 


35619 


188283 


3 


37039 


195053 


4 


44244 


215022 


5 


48459 


258595 


6 


50430 


267086 



Force 



0.3B 
0.375 

u 

« 0.37 
a 

^0.365 

o 0.35 
Eh 

0.355 
0.35 

0.345 
0.34 




'Force 



3 4 5 

Refinement Level 



Figure 4: Variation of total force with respect to refinement level 

V. Results And Discussion 
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Figure 4(a): Static pressure magnitude contour for 0° rotor angle 
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Figure 4(b): Static pressure magnitude contour for 60° rotor angle 
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Figure 4(c): Static pressure magnitude contour for 120° rotor angle 
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Figure 4(d): Static pressure magnitude contour for 1 50° rotor angle 



The helical Savonius rotor with three bladed 90 D twist angle, the contour plot analysis at 0 D , 60 D , 120 D and 150* 

rotor angle for one revolution of the rotor, the contour plot analysis gives an idea about the flow physics of a wind rotor and 
its power generation mechanism. In the present study, relative velocity magnitude (velocity of the rotor relative to wind) and 
static pressure contours are analyzed. Fig- 5(a) to Fig- 5(d) show the velocity contours of the three-bucket helical Savonius 
rotor for different rotor angle conditions of 0 , 60 , 120 L and 150 L . The velocity magnitude contour shows that the 

maximum velocity magnitudes near the chord end of advance bucket and decrease downstream side of the rotor. For 0 D 

rotor angle velocity obtained 2.98x1 0 : m/sec to 1.49m/sec, for 60 D rotor angle 2.86x1 0 1 m/sec to 1.43m/sec, for 120 D 

j p. j 

rotor angle 2.98x10 m/sec to 1.49m/sec and at 150 - rotor angle 2.95x1 Q- m/sec to 1.47m/sec at XZ iso-plane. Fig- 4 (a) to 

Fig- 4 (d) show the static pressure flow separation near to the bucket wall and Fig- 6 (a) to Fig- 6 (d) show the static pressure 
contours of the rotor for the aforesaid rotor angle conditions. These pressure contours show a decrease of static pressure from 
the upstream side to the downstream side across the rotor, which results in useful for evaluate the performance of rotor. At 
0 D rotor angle Fig- 6(a), static pressure decreases from 3.67x10* Pascal to -3.97x 10 _ Pascal from upstream side to 

p. ■ wy. 

downstream side. At 60 rotor angle Fig- 6(b), static pressure decreases from 4.36x1 0 _ Pascal to -3.56x 10 _ Pascal from 

p. n. 

upstream side to downstream side. At 120" rotor angle Fig- 6(c), static pressure decreases from 3.83x10* Pascal to -3.82x 
10 _ Pascal from advanced bucket to downstream bucket. At 150 D rotor angle Fig- 6(d), static pressure decreases from 
3.97x1 0 : Pascal to -4.004x 10 : Pascal from upstream bucket to downstream bucket at XZ iso-plane. 
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Figure 5 (a): Velocity contour of 0° rotor angle at XZ iso-plane 
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Figure 5 (b): Velocity contour of 60° rotor angle at XZ iso-plane 
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Figure 5 (c): Velocity contour of 120° rotor angle at XZ iso-plane 
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Figure 5 (d): Velocity contour of 150° rotor angle at XZ iso-plane 
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Figure 6 (a): Static pressure contour of 0° rotor angle at XZ iso -plane 
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Figure 6 (b): Static pressure contour of 60° rotor angle at XZ iso -plane 
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Figure 6 (c): Static pressure contour of 120° rotor angle at XZ iso-plane 
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Figure 6 (d): Static pressure contour of 150° rotor angle at XZ iso-plane 
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vi. Conclusion 

Flow physics analysis clearly shows that the static pressure decreases from the upstream side to the downstream 
side across the rotor, which results in useful lift for the rotor. So the maximum static pressure obtained near to the advanced 
bucket at 0 D , 60 D , 120 D and ISO* rotor angle are 367 Pascal, 436 Pascal, 383 Pascal and 397 Pascal. Again for velocity 

contour analysis the maximum velocity obtained near to the advanced bucket chord at 0 D , 60 D , 120" and 150 11 rotor angle 
are 22.3 m/sec, 25.7m/sec, 23.8m/sec and 23.6m/sec. So it can conclude, there is no sharp increase static pressure between 

Pi Pi Pi 

rotor angles of 60 to 120 L . At 60 rotor angle maximum static pressure and velocity magnitude obtained 436 Pascal and 

25.7m/sec which is marginally higher compare to other rotor angles, would ensure improved aerodynamic torque of the rotor 
as a whole during of the power stroke in one revolution of rotor. Flow physics analysis also shows the effect of positive 
static torque at variation of rotor angle. For 90° twist three buckets helical Savonius rotor at 150° rotor angle maximum 
positive static pressure distributions near the chord ends is obtained 3.97x1 0 : Pascal to 3.06x10 ~ : Pascal from upstream 
side to downstream side which gives the positive static torque. So it can be concluded from flow physics analysis that three 

Pi Pi 

bucket 90" twist angle helical Savonius rotor at 60 rotor angle maximum lift force during power stroke was evaluated and 
maximum positive static torque was evaluated at 150 c rotor angle. 
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